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Abstract 

The dynamic relaxation process for the (2+l)-dimensional SU(2) lattice 
gauge theory at critical temperature is investigated with Monte Carlo meth- 
ods. The critical initial increase of the Polyakov loop is observed. The 
dynamic exponents 9 and z as well as the static critical exponent P/u are 
determined from the power law behaviour of the Polyakov loop, the auto- 
correlation and the second moment at the early stage of the time evolution. 
The results are well consistent and universal short-time scaling behaviour of 
the dynamic system is confirmed. The values of the exponents show that the 
dynamic SU(2) lattice gauge theory is in the same dynamic universality class 
as the dynamic Ising model. 
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In recent years lattice gauge theory has continuously been developing. Numerical simu- 
lations have been extended to the gauge theory with fermions and at finite temperature. Up 
to now, however, most of simulations are devoted to the equilibrium state. Non-equilibrium 
processes have not been studied so deeply. In this letter, we report the first numerical sim- 
ulation of the short-time dynamics for the (2+l)-dimensional SU(2) lattice gauge theory at 
the critical temperature. 

Recently great progress has been achieved in critical dynamics for spin systems. For long 
it was believed that no universal behaviour would be present in the short-time regime of 
critical dynamics. However, for the critical relaxation process starting from an initial state 
with very high temperature and small magnetization, it was recently argued by Janssen, 
Schaub and Schmittmann ^ with renormalization group methods that there exist univer- 
sality and scaling even at macroscopic early times, which sets in right after a microscopic 
time scale tmic- Based on the scaling relation it was predicted that at the beginning of the 
time evolution the magnetization surprisingly undergoes a critical initial increase 

M{t)r^mot^ (1) 

where ^ is a new dynamic exponent. 

Numerical simulations support the above predictions. The critical initial increase of the 
magnetization in Eq. (|l]) was observed for the Ising model and the Potts model and the 
exponent 6 was directly measured [0,0. The scaling relation and universality are confirmed 
IQJ^J^,^-^]. The microscopic time scale t^ic is around 5 — 60 Monte Carlo time steps, depend- 
ing on the observables and the microscopic details. If one argues that a Monte Carlo time 
step is a typical microscopic time unit, such a result is reasonable. Compared with in the 
long-time regime, the scaling variable t in short-time dynamics plays a more important role 
in understanding the scaling behaviour, especially in numerical simulations. The investiga- 
tion of the universal behaviour of the short-time dynamics not only enlarges the fundamental 
knowledge on critical phenomena but also, more interestingly, provides possible new ways 
to determine all the dynamic exponents as well as the static exponents from the short-time 
dynamics, either based on the power law behaviour of the observables at the beginning of 
the time evolution |]3|j9|,H, |T0|] , or on finite size scaling |11,12|. An appealing feature of these 



methods is that they may be free of critical slowing down 

At this stage, it is natural to ask whether such an investigation may be generalized to 
field theory. This would be very important for the understanding of the non- equilibrium 
properties of field theory and numerical simulations of lattice gauge theory. On the other 
hand, the two-dimensional Ising and 3-state Potts model are known to be the simplest 
models presenting critical phenomena. They are known to have quite clean behaviour in 
many respects. One may wonder whether nice universal short-time behaviour is special for 
these simple systems. For example, it could be that the microscopic time scale tmic for more 
complicated systems is so big such that it is comparable with the macroscopic time scale. 
Then we can not observe any universal behaviour in the macroscopic short-time regime. 

As a first approach to the dynamic gauge theory, we will numerically investigate the 
(2-|-l)-dimensional SU(2) lattice gauge theory at finite temperature with a dynamics of 



model A [|1^] even though there exists so far not any analytical investigation. The motivation 
to choose this model is that the deconfining phase transition observed in this model is the 
second order. On the other hand, it is expected that in equilibrium the (2+l)-dimensional 
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SU(2) lattice gauge theory is in the same universahty class as the two-dimensional Ising 
model Ijl^. This has also been numerically observed |15,|l6|. It is interesting to investigate 



whether the dynamic SU(2) lattice gauge theory is also in the same dynamic universality 
class as the dynamic Ising model. 

The (2+l)-dimensional SU(2) lattice gauge theory is described by the Hamiltonian 

H=-^,T.UP (2) 
y p 

where the index P indicates the sum over all the fundamental plaquettes 

Up = U^{x)U^{x + fi)Ul{x + fi + u)Ul{x + z>) (3) 

with /i, V denoting the directions in the (2+1) dimensional space. As a finite temperature 
theory, the lattice size should be taken to be N'^xNq. Here A^^o corresponds to the inverse 
temperature. Since the theory is super-renormalizable, the continuum scaling law is of the 
simple form NqT ~ g~'^. For this SU(2) lattice gauge theory in equilibrium there exist 
already rather good numerical results. For a lattice with = 64 and A"o = 2, Christensen 
and Damgaard obtained the critical point A/ g1 = 3.39 and the exponent (3 = 0.120(8) ||15|| , 
while Teper got the critical point A/g"^ = 3.47 and the exponent u = 0.98(4) |jT6|. Compared 
with the exact values f3 = 0.125 and u = 1.0 for the two-dimensional Ising model, the 
numerical results of the critical exponents for the SU(2) lattice gauge theory support that 
both models are in the same universality class. Values for the critical point from the two 
groups of authors show some small difference. The determination of the critical point to 
a very rigorous level is at the present stage still difficult. For simplicity, in this paper we 
will take the simple average 4/(7^ = 3.43 of the above given two values as an input. Our 
numerical results show that this value for the critical point is very close to the real one. 

In order to simulate a critical relaxation, one should first prepare an initial state. The 
magnetization in the SU(2) lattice gauge theory is defined as the globally averaged Polyakov 
loop 

M{t) = ^T.<W^it)>, (4) 

i 

where Wi is the Polyakov loop at site i which locates in the two-dimensional lattice with 
lattice size A^^. The average < ... > is over the random forces and the independent initial 
configurations. Compared with the Ising model the Polyakov loop Wi plays the role of the 
Ising spin Si [|. Following the idea of Janssen, Schaub and Schmittmann [|l|, the initial 
state should have zero spatial correlation length and small initial magnetization. To get 
initial configurations with zero spatial correlation, we just remember that for fixed A'o, the 
center symmetry Z2 of the SU(2) theory is broken for large 'i/g'^, while it remains unbroken 



^However, Wi in the SU(2) theory is a real variable defined in the interval [—1, 1] while the spin Si 
takes integers ±1. Therefore the SU(2) lattice gauge theory should contain more physical contents. 
But in this paper we will not discuss this in details. 
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for small 4/(7^. Therefore, for the initial state we take the coupling 4/(7^ = 0. A non- 
zero initial magnetization can be achieved in many ways. A natural one is to introduce an 
initial external magnetic field h. Summarizing, the initial configurations can be generated 
by the initial Hamiltonian Hq = hJ^iWi. Final procedure is to adjust the configuration so 
generated as to give the initial magnetization sharply. Different methods can be found for 
this sharp preparation in literature in statistical physics [|],^^,||,|1^,|18|] . 

After an initial configuration is prepared, the system is suddenly quenched to the critical 
temperature with the Hamiltonian in Eq. (^) and then released to evolve with a dynamics 
of model A. In this paper we adopt the heat-bath algorithm for the dynamic evolution. 
We stop the update of the dynamic system at a reasonable time, which is typically some 
hundreds Monte Carlo time steps, and repeat the process. The average is taken over different 
independent initial configurations and random numbers. The total sample for the lattice 
size up to = 64 is 9 600 while for A^ = 128 is 8 000. Errors are estimated by dividing the 
total sample into three or five groups. 

In Fig. |I|, time evolution of the magnetization is plotted in double-log scale for different 
initial magnetization mo and lattice sizes A^. The solid line above is the magnetization profile 
for mo = 0.04 and A^ = 64, while the dotted and dashed line are those for A^ = 16 and 32 
respectively. The solid line below is the time dependent magnetization with mo = 0.02 
and A^ = 64. From the figure one can realize that in the first Monte Carlo time step the 
magnetization drops. For example, for M{t = 0) = mo = 0.04 with lattice size A^ = 64 
it drops to M(t = 1) = 0.026. Such a dropping within the microscopic time scale tmic 
is a typical non-universal behaviour which essentially depends on the microscopic details. 
Similar phenomena have also been observed for the two-dimensional Potts model and the 
XY model 0,^. With the Metropolis algorithm for the Potts model, the magnetization 
decreases continuously even up to around 10 Monte Carlo time steps. From Fig. |T], we 
clearly see that after one Monte Carlo time step, the magnetization indeed increases and for 
big enough lattice sizes it is quickly stabilized to the universal power law behaviour given 
in Eq. (|l]). The microscopic time scale is tmic ~ 20. For mo = 0.04, the magnetization 
profiles for A^ = 32 and 64 do not show a big difference. The finite size effect for A^ = 64 
is already quite small. From the slope of the curves in Fig. |1|, one can measure the critical 
exponent 6. For mo = 0.04 and A^ = 64, from a time interval [20, 250] we obtain the exponent 
6 = 0.192(2). This value is very well consistent with that for the two dimensional Ising model 
. The comparison can be done in Table |I|. In Table |, the value of the exponent f3/i> for 
the Ising model is exact while those of the exponents 6 and z are taken from the literature 
0. Rigorously speaking, the critical exponent 6 is defined in the limit mo = 0. A numerical 
measurement in this limit is practically not possible. Therefore, a linear extrapolation to 
the limit mo = from finite mo should in principle be carried out . For this reason, we 
have also performed the simulation for mo = 0.02 and A^ = 64. Since mo now is smaller, 
the fluctuations become bigger. The measured exponent is 6' = 0.186(12). Within the errors 
we can not distinguish the results for mo = 0.02 and mo = 0.04. Therefore in this paper a 
linear extrapolation will not be performed. 

Now we set mo = and proceed to measure the auto-correlation A(t) and the second 
moment M^'^\t). From a careful scaling analysis for Ising-like systems |T9|,p|, we expect the 
short-time dynamic scaling behaviours 
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A{t) = ^ E < w,m^{^) >- 1'-'/' (5) 



M(2)(t) = ^ < (E ^^it))' t^'-'^^'^/' (6) 

i 

Important is here that the exponent 6 also enters the auto-correlation A{t). Actually the first 
numerical estimate of the exponent 6 for the Ising model is from the measurement of A{t) 
by taking z as an input 0,^. In this case, however, the error induced by z is relatively big 
since the dynamic exponent z is usually not known so rigorously and the ratio d/z is much 
bigger than 6. In contrast to this, with the exponent 6 obtained from the initial increase of 
the magnetization, one can get a rather accurate value for the exponent z from the auto- 
correlation. Compared with the traditional measurement of z from the exponential decay 
of the auto-correlation in the long-time regime, to some extent, our short-time dynamic 
approach is free of critical slowing down. Due to the large time correlation length, for big 
lattice sizes the traditional measurement is very difficult since independent configurations 
can hardly be generated. However, in our short-time dynamic approach the measurement is 
always carried out in the short-time regime, where the system rather rapidly converges to 
the universal power law behaviour as the lattice size increases. In principle, we do not have 
the problem how to generate independent configurations efficiently. On the other hand, in 
the dynamic approach the average is really the sample average rather than the time average 
based on the ergodicity assumption. This fact may show its merit in future. 

Here we should stress that the dynamic exponent z is originally defined in the long-time 
regime of the dynamic evolution. It is a conceptual progress that we can determine it from 
the universal short-time behaviour. More importantly, not only the dynamic exponent z but 
actually also the static exponents, which are defined in equilibrium state, can numerically 
be measured from the short-time dynamics, e.g. from the second moment M(^^(t). 

To measure the auto-correlation A{t) and the second moment M*^^)(t), we have performed 
the simulation of mo = with a lattice size = 128. In Fig. |, A{t) and M^'^\t) are plotted 
in double-log scale with the solid line and dotted line respectively. After a microscopic time 
scale tmic ~ 50, nice power law behaviour is observed. In the numerical simulations for the 
Ising model and Potts model, it is known that the microscopic time scale for the second 
moment is somehow longer than that for the magnetization |^. This is also the case for the 
SU(2) lattice gauge theory. We have performed a power law fit in the time interval [60, 250] 
and obtain the exponent 9-d/z = -0.745(12) while {d-2p/u)/z = 0.825(14). These results 
are consistent with those of the two-dimensional Ising model, 6 — d/z = —0.737(01) and 
{d — 2P/h')/z = 0.817(07). Taking the exponent 9 as an input, we can calculate the values for 
the dynamic exponent z and the static exponent /5/z/. All these results and, for comparison 
those of the Ising model are given in Table |. Even though the dynamic exponent z has been 
known for a long time, for the lattice gauge theory our measurement is the first reliable one. 
The measured value jS/u = 0.120(18) is also in agreement with (3/v = 0.120(08) obtained in 



reference |T5|, and with the exact value iijv = 0.125 for the two-dimensional Ising model. 

In conclusion, we have investigated the universal short-time behaviour of the (2+1)- 
dimensional dynamic SU(2) lattice gauge theory. The critical initial increase of the Polyakov 
loop is observed. The dynamic exponents Q and z as well as the static critical exponent fijv 
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are determined from the power law behaviour of different observables at the early stage of the 
time evolution. The values are consistent within the errors with those of the two-dimensional 
Ising model. The results strongly support that there exists an universal short-time scaling 
behaviour for the dynamic SU(2) lattice gauge theory, and also suggest that the (2-1-1)- 
dimensional dynamic SU(2) lattice gauge theory and the two-dimensional Ising model are 
in the same dynamic universality class. 
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TABLES 





e 


z 


13/v 


SU(2) 


0.192(02) 


2.135(27) 


0.120(18) 


Ising 


0.191(01) 


2.155(03) 


0.125 



TABLE I. The exponents 6, z and 2/3/i/ measured from the short-time dynamics with the 
heat-bath algorithm. 
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FIGURES 
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FIG. 1. The time evolution of the magnetization for different lattice sizes N and initial 
magnetization mo with the heat-bath algorithm is plotted in double-log scale. 
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FIG. 2. The auto-correlation A{t) and the second moment M^'^^t) for the lattice sizes = 128 
and initial magnetization mo = obtained with the heat-bath algorithm is plotted in double-log 
scale. 
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